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Resonances in the magnetic decoupling curves for the spin 
relcixation of dense alkali-metal vapors prove that much of the 
relaxation is due to the spin-axis interaction in triplet dimers. 
Initial estimates of the spin-axis coupling coefficients for the 
dimers (likely accurate to a factor of 2) are \X\/h = 290 MHz 
for Rb; 2500 MHz for Cs. 

Hot, dense alkali- metal vapors are used to polarize the 
nuclear spins of noble gases, especially ^He and ^^^Xe, 
by spin exchange optical pumping |l|. Applications of 
these "hyperpolarized" noble gases in medical imaging 
[^ and fundamental physics ||^ are of considerable cur- 
rent interest. Nevertheless, the basic physical limits to 
the efficiency of the spin-exchange process are still not 
fully understood. Under typical operating conditions of 
Rb-^He spin exchange polarizers (saturated Rb vapor at 
200° C in 8 amagat of ^He) about 45% of the Rb spin 
depolarization is due to Rb-Rb collisions B . Here we re- 
port conclusive experimental evidence that much of the 
spin relaxation in dense alkali-metal vapors comes from 
the spin- axis interaction in triplet dimers. Our measured 
spin-axis coupling strengths also impact the field of ul- 
tracold collisions, where the spin-axis interaction causes 
trap- loss and affects the widths and positions of Feshbach 
resonances ^g- 

In 1960, Anderson et al. H] showed that the dominant 
effect of binary collisions between alkali-metal atoms is 
exchange of the electron spins. In 1974 Gupta et al. ||] 
showed that the spin angular momentum of the alkali- 
metal atoms was freely exchanged with the nuclear spin 
of the singlet dimers, for example, Rb2 and Cs2 (where 
nuclear quadrupole interactions cause some spin depo- 
larization). In 1980, Bhaskar et al. H showed that Cs- 
Cs interactions destroy spin at about 1% of of the spin- 
exchange rate. The rate coefficient for Cs-induced re- 
laxation showed little dependence on buffer gas pressure 
or species. Relaxation mechanisms in singlet or triplet 
dimers were expected to depend strongly on pressure. 
Consequently, Bhaskar et al. ||] proposed that the relax- 
ation was due to binary collisions between alkali-metal 
atoms. They pointed out that although there should be 



a spin-rotation interaction, 14r — 7N • S in the triplet 
dimers - similar to that postulated by Bernheim |10[| - 
the coupling coefficient 7 was probably too small to ac- 
count for relaxation in binary collisions. Instead, they 
suggested that the source of the relaxation was a spin- 
axis interaction of the classic form mU 
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Here S = Si + S2 is the total electron spin of the col- 
liding pair, and ^ is a unit vector along the direction of 
the internuclear axis. The spin-axis coupling coefficient 
A — X{R) arises from both the interaction energy of the 
electrons' magnetic dipole moments and the spin-orbit 
interaction in second order |12| . 

Because Knize's [|l^ measurements also showed little 
dependence of the Rb-induced relaxation on He or N2 
pressure, it came as a surprise when in 1998, Kadlecek, 
Anderson and Walker Mm discovered that the relaxation 
could be substantially suppressed by externally applied 
magnetic fields of a few thousand Gauss. If the interac- 
tion were due to pairs of alkali-metal atoms interacting 
on the triplet potential curve, this field dependence sets 
a lower limit on the correlation time of 40 ps, the in- 
verse of the electron Larmor frequency at the decoupling 
field and a time much longer than the duration of a bi- 
nary collision. Similar magnetic decoupling curves for K 
and Cs vapors were soon observed in our laboratories at 
Wisconsin and Princeton. 

Although we still do not understand the relaxation at 
He or N2 pressures of an atmosphere or more, recent low- 
pressure experiments leave no doubt that an important 
part of the spin relaxation in dense alkali-metal vapors 
comes from the spin-axis interaction (1), acting in triplet 
dimer molecules - even though the triplet dimer density 
is no more than 10~^ that of the monomers. The key 
experimental observation is the existence of resonances 
in the magnetic decoupling curves, resonances which are 
predicted from the spin-axis interaction (1), and which 
cannot be produced by the spin-rotation interaction or 
by spin- 1/2 species such as trimer molecules. 

A representative arrangement for measuring relaxation 
transients in Rb is shown at the top of Fig. 1. The alkali 
vapor and N2 buffer gas were contained in 1 inch spherical 
glass bulbs that were heated in an oven placed between 
the pole faces of a 10 kC electromagnet. We measured 



the number density of Rb atoms using ^i^^^y rotation 
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Faraday rotation of linearly polarized light both to mea- 
sure the Cs density and to monitor the spin-polarization. 
Representative measurements of 7 in *^Rb are shown 
in Fig. 2(a), where one can recognize 3 well- resolved res- 
onances on the magnetic decoupling curve at magnetic 
fields of approximately 900, 1500 and 2100 Gauss. Sim- 
ilar representative data for ^^'^Cs are shown in Fig. 3, 
where the theoretical curve has 7 poorly resolved reso- 
nances, which nevertheless provide a good fit to the ex- 
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Resonances due to SPIN-AXIS Interaction 
X/h = 290 MHz 
\1/x = 3.36 X 109 s-1 



FIG. 1. Top. Experimental apparatus for measuring the 
spin relaxtion of Rb vapor. Bottom. Representative transient 
decay curve obtained with the apparatus; at late-times the 
transient decay is characterized by a single exponential time 
constant. 

An intense pulse of 7947 A circularly polarized pump 
laser radiation built up the spin polarization for a short 
(10 ms) interval. The pump laser was then blocked by 
the chopper wheel, and photodiodes monitored the trans- 
mission of a weak linearly polarized probe laser through 
the cell. The probe beam was resolved into its two 
circular components, giving positive and negative helic- 
ity signals V±(t) for sampling times t. The signal ra- 
tio I{t) — V+{t)/V-{t) greatly suppressed probe-laser 
noise. The difference in the attenuation coefficients of 
the two circular components of the probe laser is pro- 
portional to the spin polarization P (for P -^ 1) of the 
vapor, so lnl{t) = a + bP{t), where a = ln/(c!o) and 
b are constants. Fig. 1 shows a typical decay transient 
lnln7(t)//(oo). For early times, the decay is not expo- 
nential since the faster eigenmodes m^ of spin relaxation 
are still contributing to the polarization P. However, 
at late times there is a time interval, several e-foldings 
in length, over which lnln/(i)//(cxD) is nearly a straight 
line with a slope —7. Thus, when the polarization has 
decayed for a sufficiently long time, it is characterized by 
a single, exponential late-time decay rate 7. 

The late-time decay rates of polarized Cs vapor in N2 
gas were measured as described previously llj], using 
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FIG. 2. Top. Late-time decay rate of ^'^Rb versus mag- 
netic field B. Resonant enhancements of the decay rate 
are observed at the fields Bi = 915 G, B2 = 1526 G 
and B3 = 2136 G, in agreement with the predicted values 
Bi = A{21 + l)/4(7s/is. The cell contained 0.061 amagat of 
N2 gas. The cell temperature was 220 C. The solid line is 
a theoretical fit assuming a field-independent contribution to 
7 of 97 sec~^ from binary collisions and diffusion to the cell 
walls, a field-dependent contribution from quadrupole interac- 
tions in singlet dimers, calculated as outlined in Ref. Iql, and 
a contribution from a spin-axis interaction in triplet dimers, 
calculated according to Eq. (kl). Bottom. Energy levels for a 
representative spin space with 1 = 2. The level anticrossings 
are responsible for the resonances in the magnetic decoupling 
curves. 

Also shown in Fig. 2(b) are representative energies of 
the spin sublevels of the triplet dimer molecules of *^Rb 
for a value of the total nuclear spin quantum number 
1 = 2. The resonances in the relaxation rate occur at the 



same magnetic fields as anticrossings in the energy level 
diagrams. The anticrossings are the result of the spin- 
axis interaction (1), acting in a total spin Hamiltonian 
for homonuclear triplet dimcrs of the form 

H=^S-I + gsf^BBS, - ^S . (3nn - 1) • S. (2) 

Here A is the magnetic-dipole hyperfinc coupling coeffi- 
cient for the free alkali atoms {A{^'^Rh)/h = 3, 417.4 MHz 
and A{^^^Cs)/h = 2,298.2 MHz). The operator for the 
total nuclear spin, I = Ii + 12 is the sum of the nuclear 
spin operators Ii and I2 for the two atoms. The operator 
I • I commutes with (2), and for homonuclear dimers its 
eigenvalues /(/ + 1) can have / = 0, 1, 2, . . . , 2I„i, where 
Im ^ Ii = I2 is the nuclear spin quantum number of a 
monomer. The Hamiltonian (2) can therefore be diago- 
nalized in a spin subspace / of dimension 3(2/-(- 1) to find 
eigenfunctions \i) = \Ii) and eigenvalues Ei = En- For 
A = there is a curious degeneracy in the spin subspace 
/ ^ 0, with 2/ -|- 1 sublevels crossing at the magnetic 
field Bi = A{21 + l)/4gsfJ.B- The last term in @ is the 
spin-axis interaction (1), where we assume that the rota- 
tional angular momentum N of the dimers is large enough 
that we may treat it as a classical vector, pointing along 
the unit vector n = N/N, and make the replacement 
(3CC-l)^-(3nn-l)/2. 

We estimate the relaxation caused by the interaction 
(2) with the following simple model. Let the total longi- 
tudinal spin operator of the triplet dimer be Gz — Sz+h- 
When triplet dimers are formed in dense alkali-metal va- 
pors with low spin polarization, the spin density matrix 
of the dimers will be very nearly p(0) = e^'~^' /Z, where 
/? « 4P <IC 1 is the spin-temperature parameter 
The partition function is very nearly Z = 3(2/™ -I- 1) 
The spin of the dimer at the instant of formation is 
therefore {Gzo) = "^i^m + Im + l)/3/3, where we have 
neglected non-linear terms in the small parameter /3. 
The mean change in the density matrix by the time it 
breaks up again is Ap = {U p{0)U~^) — p(0), where the 
time evolution operator is U = g-*^*/'*^ and the angle 
brackets (• ■ •) denote an average over all directions of n, 
and also an average over the probability e~*^'^dt/T that 
the dimer will break up in the time interval between t 
and t -\- dt. The collisional breakup rate of the dimers 
is l/r. The mean change in the dimer spin is there- 
fore (AGz) = Tr{G^Ap} = -W{Gzo), where the spin- 
destruction probability is 



W = 2Y^ 



mGz\J)\'icO^. 



(2/,„ + 1)2 [(2/„, + 1)2 + 3][1 + Kt)2 



(3) 



and the Bohr frequencies are ujij = {Ei — Ej)/h. An 
average over the directions of n is to be understood in (||). 
Almost all of the spin is carried by monomers of number 
density [Ak] {e.g., [Ak] = [Rb] or [Ak] = [Cs]), which have 
a spin density per unit volume [Ak]/3(4Z2^ + 4/„j + 3)/12. 
The spin loss rate per unit volume from triplet dimers 



of number density [Ak2] is [Ak2](AG'2)/T, from which 
one obtains the relaxation equation —{d/dt)hi(3 = ^ = 
K[Ak], where the rate coefficient is 



(2/„ 



f2/„ 



1)= 



WK, 



(4) 



The chemical equilibrium coefficient of the triplet dimers 
is IC — [Ak2]/[Ak]2. We have used theoretical triplet 
dimer potential curves of Krauss and Stevens ||l8| to es- 



timate /C(Rb) = 330 A3 ai 



d /C(Cs) = 650 A3 at repre- 



sentative temperatures of 2!!y°C and 180°C respectively. 
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FIG. 3. Late-time field-dependent decay rate of ^^''Cs 
versus magnetic field and pressure at 180°C, [Cs]= 6.6 x 10^* 
cm~^. The cell temperature was 177 C. The fits to the data 
use the reorientation model with negative A. 

Non-vanishing values of A lift the degeneracy at the 
field Bi and convert the level crossings into anticrossings, 
while making the corresponding matrix element {i\Gz\j) 
non-zero, as is illustrated in Fig. 2. It is these anticross- 
ings of the energy levels that are responsible for the res- 
onances, most clearly visible in Fig. 2. Remarkably, the 
spin-rotation interaction does not break the degeneracy 
so the spin-axis interaction must be responsible for the 
relaxation resonances. Since there are 2Im nonvanish- 
ing values of / for atoms with the nuclear spin quantum 
number /„j, we expect there to be 3 resonaces for *^Rb 
with Im — 3/2, as observed, and 7 resonances, at fields 
between 600 G and 3100 G for ^^^Gs with /,„ = 7/2, 
which fits well the observed field dependence shown in 
Fig 3. 

We fit the experimental data to theoretical curves that 
include a field-independent part due to diffusion to the 
cells walls and binary collisions, a field-dependent nu- 
clear quadrupole relaxation in singlet dimcrs (negligible 
for Cs), and relaxation due to the spin-axis interaction 



in triplet dimers. We expect a small portion of the re- 
laxation in triplet dimers to arise from the spin-rotation 
interaction, and we may estimate the magnitude of this 
relaxation mechanism by assuming that the spin-rotation 
interaction of Rb-Rb triplet dimers is the same as that 
that of Rb-Kr, and that the spin-rotation interaction for 
Cs-Cs is the same as that for Rb-Xe, both of which are ex- 
perimentally known from the work of Bouchiat et al. [[l9| . 
These estimates should be fairly reliable since the spin- 
orbit interaction that produces the spin-rotation coupling 
PI varies only slightly from Rb to Kr or Cs to Xe. Using 
these estimates for 7N • S, we find that the addition of 
the spin-rotation interaction causes a few percent change 
in the predicted relaxation rates. 

In the model described above, the relaxation due to 
triplet dimers dependends on two parameters: A and 
the molecular lifetime r, assuming the theoretical val- 
ues for /C are accurate. The molecular lifetime should 
be inversely proportional to pressure, making it natu- 
ral to introduce a cross-section a — (t[N2]v)~^. Our 
Rb data is well-described over the limited pressure range 
from 50 Torr to 300 Torr by the values |A| = 290 MHz, 
and a — 290 A^. The latter is somewhat larger than 
the typical 150 A^ breakup cross sections deduced from 
magnetic decoupling studies of RbKr ||T^ . The sign of A 
cannot be determined from the data. 

Since there is little contribution to Cs relaxation from 
singlet molecules, relaxation from Cs-Cs triplet molecules 
can be observed over a wider range of pressures. Fig. 3 
shows a sampling of our spin-relaxation data for Cs as 
a function of magnetic field at a variety of N2 pressures. 
In this case the data are well-described by A = —2.9 
GHz, cr = 98 A^, /C = 350 A^, or, assuming the op- 
posite sign of A, A = 3.4 GHz, ct = 82 A^, /C = 307 
A-^. The deduced values of K. are about a factor of 2 
smaller than predicted from the potential curves (which 
are believed to be fairly reliable). We have also fit the 
data with a model that allows for coUisional reorienta- 
tion of the triplet molecules without breakup, making 
the assumption that a spin-temperature distribution de- 
scribes the dimcr density matrix after reorientation. In 
this case r in Eq. a is replaced by the coherence time Tc, 
and W = Wt/{tc + (r - Tc)W) replaces W in Eq. |. In 
this case we obtain A = -2.13 GHz, CTc = 81 A^, ct = 42 
A2, K, = 584 A3, or A = 2.71 GHz, ac = 93 A^, cr = 55 
A^, K = 431 A^. We emphasize that single values of 
A and ac accurately represent the magnetic field depen- 
dences at 9 pressures (not all shown in Fig 3), covering 
nearly two decades in pressure. 

Spin-axis coupling in Cs triplet dimers is also of cur- 
rent interest for studies of ultracold collisions g since it 
produces loss of spin-polarized atoms in magnetic traps, 
and affects the widths and positions of Feshbach reso- 
nances in Cs collisions [g. The magnitude and sign of 
A, as well as its dependence on interatomic separation R 
were predicted by Mies et al. [|2|. Previous analysis g] 
led to the conclusion that in order to agree with exper- 
iment the second-order spin-orbit calculated by Mies et 



al. needed to be scaled by a factor of Sc — 4.0 ± 0.5, 
recently updated to 3.2 ±0.5 0. We find, by averaging 
\{R) over the ro- vibrational levels of the triplet state us- 
ing the Krauss and Stevens potentials ||l^ , scaling factors 
of Sc = 12 for Rb and Sc = 3.6 for Cs are needed to 
reproduce our data. 

In conclusion, the resonances observed in the magnetic 
decoupling curves of Rb, coupled with the observed pres- 
sure and field dependences of Cs, prove that a major con- 
tribution to the spin relaxation of hot, dense alkali-metal 
vapors is caused by the spin-axis interaction in triplet 
dimers. Knowing this key piece of the physics should 
help to unravel the still-puzzling independence of the re- 
laxation rates on He and N2 pressures of an atmosphere 
or more |ll4[| . 
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